Screening and collective modes in gapped bilayer graphene 
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We study the static and dynamic screening of gapped AB-stacked bilayer graphene. Unlike 
previous works we use the full 4-band model instead of the simplified 2-band model. We find 
that there are important qualitative differences between the dielectric screening function obtained 
using the simplified 2-band model and the 4-band model. In particular, within the 4-band model, 
in the presence of a band gap, the static screening exhibits Kohn anomalies that are absent within 
the simplified 2-band model. Moreover, using the 4-band model, we examine the effect of trigonal 
warping on the screening properties of bilayer graphene. We also find that the plasmon modes have 
a qualitatively different character in the 4-band model compared to the ones obtained using the 
simplified 2-band model. 



Bilayer graphene pQ has many unique electronic prop- 
erties that make it an extremely interesting system. AB- 
stacked bilayer graphene (BLG) is formed by two Bernal 
stacked layers of graphene [2-4]. When placed on an insu- 
lating substrate the electrons in BLG form an ideal two- 
dimensional electron gas (2DEG) with a very high room 
temperature mobility, in particular when boron nitride 
is used as a substrate [5]. In pristine BLG the conduc- 
tion and valence bands touch at points, charge neutrality 
points (CNPs), at the corners of the Brillouin zone. At 
very low energy around these points the bands are ap- 
proximately parabolic. However, by applying a perpen- 
dicular electric field a band-gap (A) can be opened, and 
tuned [6J. Moreover, recent experiments [7] have shown 
strong evidence that at low temperatures and dopings 
the electrons in BLG might be in a spontaneously bro- 
ken symmetry state |H]. All these facts make BLG an 
extremely interesting system both from a fundamental 
point of view and for its possible technological applica- 
tions. As a consequence the accurate knowledge of the 
electronic properties of BLG is of great interest. 

One of the most important physical quantities for char- 
acterizing the electronic properties of a system is the di- 
electric function e(q, uo). This quantity determines the 
effective, screened, Coulomb interaction among the elec- 
trons in the system and is therefore essential for the 
calculation of all the electronic properties. There is 
strong evidence that in most BLG samples charge im- 
purities close to the surface of the substrate, or placed 
between the substrate and the BLG layer, are the domi- 
nant source of scattering [4] . In this situation the knowl- 
edge of the static dielectric function, e(q, uo = 0), is es- 
sential to calculate the dc conductivity. Moreover, in 
the case of magnetic adatoms placed on BLG, the static 
polarizability determines the effective Ruderman-Kittel- 
Kasuya-Yoshida (RKKY) interaction between the mag- 
netic adatoms. The dynamic dielectric function deter- 
mines the optical properties of the system and the col- 
lective electronic modes, plasmons. It is therefore evident 
that the knowledge of the correct form of e(q, u) is nec- 
essary to characterize the electronic properties of BLG. 



Previous works [9] [10] have studied the case of gapless 
BLG (and gapless single and multilayer systems [H]). I n 
the presence of a gap some of the symmetries that sim- 
plify the calculation of the response functions in gapless 
BLG disappear. In part for this reason the only results 
available [12 for the dielectric function in gapped BLG 
were obtained using a simplified effective low energy 2- 
band model p~3j [14] . This model neglects features of the 
band structure of BLG that can strongly affect the re- 
sponse function, especially when A In particular, it 
neglects the fact that in the presence of a band-gap the 
bands, at low-energy, acquire a characteristic "sombrero" 
shape [13] , see Fig. [l](a). To describe these effects it is 
necessary to use a refined 2-band model [l3j[T4] or the full 
4-band model. In this work we obtain e(q, ui) for gapped 
AB-stacked bilayer graphene using the full 4-band model 
and the random phase approximation (RPA). Some of the 
qualitative differences for e(q, cj) between the full 4-band 
model and the simplified 2-band model can be recovered 
using the refined 2-band model. However, there are fea- 
tures of e(q, uj) (especially at large /c, uo and n) obtained 
using the full 4-band model that are qualitatively dif- 
ferent from e(q,cj) obtained using either the simplified 
or the refined 2-band model. In the remainder, unless 
specified, by 2-band model we refer to the simplified one. 

The 4-band continuum model Hamiltonian for BLG 
is H = -E k ^(k)^k where tf* (tf k ) is the 
4-component creation (annihilation) operator Sfr^ = 
( a k,i' & k,i' a k,2' & k,2) (^k = (ak,i,&k,i,ak,2,&k,2)) with 
&k i (&k,i), b\ ii (6k,i) the creation (annihilation) operator 
for an electron with wavevector k in layer i on sublattice 
A and B respectively, and h{k) is the matrix 

A 71 
ft(k) = —r z + hv F (k x a x + ky(jy) - — (cf x t x + a y r y ) 
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In Eq. Q a's, r's are 2x2 Pauli matrices representing the 
sublattice and layer degrees freedom respectively, vp is 
the Fermi velocity at the Dirac point of a single graphene 
layer, 71, 73 are the interlayer hopping parameters [15], 
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a = 1.42 A is the in plane carbon-carbon distance, and 
A is the band gap at k = 0. 73 ^ induces trigonal 
warping. For concreteness, we assume v F = 10 6 m/s, 
71 = 0.35eV, and 73 = (3/4)71 = 0.26eV, however the 
main features of our results do not depend on the precise 
values chosen for these parameters. 

The Coulomb interactions are described by the Hamil- 
tonian H % = (1/2 A) £ q [V+(g)/3 q p_ q + V^)d q <L q ], 

where A is the sample area, p q (<i q ) the operator for 
the sum (difference) of the densities in the two lay- 
ers, V±(q) = (V s (q) ± V D (q))/2 with V s (q) = 2^e 2 /(eq) 
the Coulomb interaction between electrons in the same 
layer and Vd = 2ne 2 (e~ qd ) / (eq) the Coulomb interac- 
tion between electrons in different layers, d = 3. 35 A the 
distance between the two layers, and e the background 
dielectric constant. We assume a = e 2 /ehv F = 0.5 and 
temperature T = 0. As long as q <C 1/d the dielectric 
function that enters the calculation of most of the elec- 
tronic quantities is the one associated with the sum of the 
densities in the two layers, e(q, uj) = e(q, oS) pp . Within 
the random phase approximation 



e(q,^) pp = 1 - V+(q)II(q,a;) 



where 

n(q,cj) pp - g 



^ r dk 



pp 



^A,k - ^A',k+q 



(2) 



) 2 fiw + e Ajk - ev >k+q + ir) 



x |(£/ k £/ k+q )A,A'| 2 



(3) 



is the polarizability (in the remainder of this paper the 
subscript pp will be understood). In Eq. Q g = g s g v = 4 
is the total spin (g s ) and valley (g v ) degeneracy, A, A' are 
the band indices, (e^,k) is the Fermi-Dirac distri- 

bution function (energy) for a particle in band A with 
wavevector k, and is the unitary matrix that diago- 
nalizes the Hamiltonian Ho. From Fig.[l](b) we see that 
the intraband wave-function overlap |(/7j[/7k+ q )A,A| 2 f° r 
the 4-band model is quite different from the one for the 
2-band model, especially when Below we present 

our results, obtained using Eq. ([2|-([3| and evaluating the 
integral on the r.h.s of Eq ^ numerically using an adap- 
tive integration routine [T6] . 

In Fig. [I] (c) the results for the static polarizability 
II(q, uj = 0) are shown for fixed doping n = 10 12 cm -2 
and different values of A. We see that for A ^ 
the results obtained with the 4-band model are very 
different from the ones obtained with the 2-band model. 
In the 2-band model II (q, 0) exhibits a Kohn anomaly 
only for q = 2k f (k F being the Fermi wavevector), 
whereas in the 4-band model (and the refined 2-band 
model) II(q, 0) exhibits Kohn anomalies also for values 
of q < 2kp. This is due to the fact that in the 4-band 
model, at low energies, the lowest bands, for A ^ 0, 
acquire a typical nonmonotonic sombrero shape. As a 
consequence in the 4-band model, for A 7^ 0, for fixed n 
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FIG. 1. (Color online) (a) Lowest conduction band for 
A = 71/2, 7i/3, 0. The solid (dash-dot) curves are obtained 
using the 4-band (2-band) model. The horizontal dashed lines 
indicate the 4-band Fermi energy for doping n — 10 12 cm -2 for 
A = 71/2, 7i/3, from top to bottom, (b) Chirality factors, 
\(UlUv) X \\ 2 , evaluated for |k| = |k ; | = k F for the 2-band 
model at A = 0, 7i/2, denoted by the black solid and dashed 
lines respectively, and the 4-band model for A = 71/2, 0. 
6 is the angle between k and k ; . For the case A = 71/2 
in the 4-band model there are three possible intraband over- 
lap scenarios we can consider: (i) k and k ; both lie on the 
Fermi surface at wavevector &f+, (ii) k and k 7 both lie on 
the Fermi surface at wavevector kF-, (in) k lies on the Fermi 
surface at wavevector while k ; lies on the fermi surface 
at kF-- (c) II(q, 0) for n = 10 12 cm -2 without trigonal warp- 
ing. Solid (dashed) curves are the results obtained using the 
4-band (2-band) model, (d) Contour plot of polarizability, 
n(g, n, uj = 0), as a function of q and doping n for A = I/271. 



(A) when A > A c = Hv F ^/7m (|n| < n c = /S 2 / {nh 2 v 2 F )) 
the Fermi surface is multiply connected. Neglect- 
ing trigonal warping for n < n c the Fermi surface 
is formed by two circumferences, of radii kp± = 

(l/hv F )yf% + A 2 /4 ± ^(7? + A2 ) - res P ec " 
tively, with e F = (l/2)^(h 4 v 4 F 7T 2 n 2 + A 2 7 2 )/( 7 2 + A 2 ) 
(see inset of Fig. [I] (b)). In this situation we can 
expect additional Kohn anomalies for values of q joining 
points on the same connected part of the Fermi surface 
and on disconnected parts of the Fermi surface. For 
71 = 350 meV and n = 10 12 cm -2 we have that 
A c « 7i/3. For A = A c the Fermi energy just touches 
the top of the sombrero. In this case we only have 
one additional Kohn anomaly for q = k F in addition 
to the q = 2kp one, see Fig. [I] (c). For A > A c the 
Fermi energy cuts the sombrero region and so we have 
Kohn anomalies for q = kp+ — fcp_, and q = 2k F - in 
addition to the one for q = 2k f+ as shown in Fig. [l] (c). 
One might expect to observe an anomaly also for 
q = k F + + k F -, however the points on the Fermi surface 
connected by this value of q have Fermi velocities with 
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the same sign and therefore the anomaly is suppressed. 
Fig. [I] (d) shows the dependence of II(q, 0) on q and n 
for A = 71/2. From this figure we see the evolution of 
the Kohn anomalies with doping, in particular we can 
observe the merging of some of the anomalies for specific 
values of the doping. 

We now consider the effects on II (q, 0) of trigonal warp- 
ing. As shown by the left panels of Fig. |2j in the presence 
of trigonal warping the energy-bands become anisotropic 
P~3j[14]. In particular, at low energies the lowest bands 
exhibit 4 degenerate minima. The modifications of the 
fermionic energy bands due to the trigonal warping are 
reflected in the polarizability, as shown by the right pan- 
els of Fig. [5] II(q, 0) becomes strongly anisotropic, the 
number and position of the Kohn anomalies becomes de- 
pendent on the direction of q. 



a) E ( Y .) b) n(q,0) 
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FIG. 2. (Color online). (Left Column) Equipotential lines 
for the lowest energy band within the 4-band model with 73 = 
37i/4 and A = 0, I/371, I/271 from top to bottom. (Right 
Column) II(q, 0) for n = 10 12 cm -2 , trigonal warping 73 = 
37i/4, and A = 0, I/371, I/271 from the top panel to bottom 
one. /cf,7 3 =o is in the limit 73 = 0. 



see that as A crosses A c the dispersion of the plasmon 
mode outside the particle-hole continuum is not modi- 
fied qualitatively. The plasmon mode inside the particle 
hole-continuum on the other hand is qualitatively very 
different for A < A c and A > A c , an effect that is not 
captured by the 2-band model [T2] . 




FIG. 3. (Color online). The left (right) column shows the 
real (imaginary) part of £RPA(q, uj) for A = 0, A = 7i/3, 
A = 71/2 from top to bottom. The plasmon dispersion is 
denoted by white curves. The boundaries for the intraband 
(interband) particle-hole continuum are indicated with black 
solid (dashed) curves. 



In the presence of trigonal warping e(q, uS) becomes 
strongly anisotropic and this is particularly evident when 
the Fermi energy cuts the sombrero region. Fig. [4] shows 
the results for e(q, uo) for different directions of q obtained 
taking into account trigonal warping. From the figure the 
strong anisotropy of e(q, ui) when 73 =^ is evident. In 
particular, we see that the plasmon mode inside the p-h 
continuum is very different for different directions of q. 

For the case with no trigonal warping in the long- 
wavelength limit q <C uo/vp using the 4-band model for 
the polarizability, up to order q 2 , we have: 



The dynamic dielectric function e(q, lj) for fixed dop- 
ing n = 10 12 cm" 2 and A < A c , A = A c , A > A c for the 
case in which 73 = (no trigonal warping) is shown in 
Fig. [3] The white lines show the plasmon dispersion, 
the black solid (dashed) lines show the boundaries of 
the intraband (interband) particle-hole continuum. We 



n(q,cj) 



<9e k 
dk 



<9e k 
dk 



(4) 



We notice that in Eq. Q there is a term proportional to 
kp- that is absent in the 2-band model. Replacing this 
expression in the equation for the RPA e(q, uS) we find 
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FIG. 4. (Color online). The left (right) column shows the 
real (imaginary) part of £RPA(q, uj) with A = 71 /2 for = 
0° , 15° , 30° . The white curves denote the plasmon dispersion, 
the black curves denote the boundaries of the particle-hole 
continuum. kF,-y 3 =o is in the limit 73 = 0. 



model does not capture the nonmonotonic band struc- 
ture, i.e. the fact that in the 2-band model in Eq. ([5| 
there is no term kp-de^/ d~k\k F _ . For n > n c this is 
due to the fact that in the 4-band model the dispersion 
is much closer to linear than parabolic as assumed in 
the 2-band model, in analogy to what happens in the 
gapless case [10]. As a consequence for very large n we 
have UJ4- band /^2-band °c rT 1 ^. This is summarized in 
Fig. [5] (b) which shows the ratio co>4_ band /^2-band be- 
tween the plasmon frequency obtained within the 4-band 
and the 2-band model as a function of n for different val- 
ues of A. Notice that in the long- wavelength limit this 
ratio (see Eq. (|6|) is independent of q and is a function 
only of n and A. 




FIG. 5. (Color online) (a) Plasmon dispersion for A = 71/2 
obtained using the full 4-band model, and the 2-band model 
for n — 2.7 x 10 12 cm -2 . (b) Ratio CJ4- band /^2-band as a 
function of doping for different value of A. For A / and 
n —> the ratio 

^4 — band /^2-band diverges. 



the plasmon dispersion: 



9 k 

-nvpaq 



<9e k 
<9k 



<9e k 
<9k 



1 1/2 



(5) 

This dispersion is very general and is valid both for n < 
n c and n > n c , in the latter case kp- =0. From Eq. ([5| 
using the appropriate expressions for kp+, kp- and ek 
we find 



^) = J^F(n,A) 

V e Ft 



(6) 



where n = h v F 7rn/ji, A = A/71, and, for the 2-band 
model, cp = 7i[n 2 + A 2 ^] 1 / 2 , F(n, A) = n, whereas 
for the 4-band model in the sombrero region, ep = 
7l [(n 2 + A 2 )y(l + A 2 )] 1 /2/ 2 , F(h, A) = [n(A 4 + 2A 2 - 
n 2 )/(A 4 + 2A 2 - n 2 + l)] 1/2 /2, and outside the sombrero 
region, e F = 71 [2 + A 2 + 4n - 2[1 + 4n(l + A 2 )] 1 / 2 ] 1 / 2 ^, 
F(n, A) = [(n/2)[l-(l + A 2 )]/[l + 4n(l + A 2 )] 1 / 2 ] 1 / 2 . In 
Fig. [5] (a) we compare the results for the plasmon disper- 
sion obtained numerically using the 4-band model with 
the ones given by Eq. (|6| for the gapped case A = 71/2 
for a given value of n. We see that the 2-band results 
differ substantially from the 4-band results. At low den- 
sities (n < n c ) this is due to the fact that the 2-band 



In conclusion, we have obtained the static and dielec- 
tric screening of gapped BLG using the full 4-band model. 
We find that the static screening obtained using the 4- 
band model is qualitatively different from the one ob- 
tained from the 2-band model. In particular in the 4- 
band model, when the gap is nonzero, the static polariz- 
ability exhibits Kohn anomalies that are not present in 
the simplified 2-band model. For the dynamic screen- 
ing we have found that the plasmon frequency within 
the 4-band model is substantially different from the one 
obtained within the 2-band model especially at low den- 
sities when We have also characterized the strong 
anisotropic properties of the static and dynamic screen- 
ing due to the trigonal warping. We find that in the 
presence of trigonal warping in gapped BLG the number 
of Kohn anomalies depends not only on the doping and 
the band-gap but also on the direction of the momentum. 
Our results, in particular the identification of additional 
Kohn anomalies, and the strong anisotropic nature of 
the screening in the presence of trigonal warping, have 
important implications for understanding of the phonon 
spectrum and the nature of the RKKY interaction in 
gapped BLG, and are therefore expected to have clear 
experimental signatures. Moreover our results apply also 
to the case in which a gap opens due to the realization of 
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a spontaneously broken symmetry state and could then 
be used to identify and characterize such a state. 
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